NUSC  Technical  Report  8689 
18  May  1990 

DTIG  FILE  COPY 

AD-A224  906 

Evaluation  of  Integrals  and  Sums 
Involving  [sin(Mx)/sin(x)]'^ 


Albert  H.  Nuttail 
Surface  ASW  Directorate 


Naval  Underwater  Systems  Center 

Newport,  Rhode  Island  /  New  London,  Connecticut 


Approved  for  public  release;  distribution  is  unlimited. 


Preface 


This  research  was  conducted  under  NUSC  Project  No.  A75215, 
Subproject  No.  ROONOOO,  "Determination  of  Concentrated  Energy 
Distribution  Functions  in  the  Time-Frequency  Plane",  Principal 
Investigator,  Dr.  Albert  H.  Nuttall  (Code  304).  This  technical 
report  was  prepared  with  funds  provided  by  the  NUSC  In-House 
Independent  Research  and  Independent  Exploratory  Development 
Program,  sponsored  by  the  Office  of  the  Chief  of  Naval  Research. 
Also,  this  report  was  conducted  under  NUSC  Project  No.  D19840, 
Subproject  No.  S0223,  "Large  Area  Array",  Principal  Investigator 
Dr.  H.  H.  Schloemer  and  Dr.  S.  H.  Ko  (Code  2133),  and  sponsored 
by  R.  Dosti,  NAVSEA-63D4. 

The  technical  reviewer  for  this  report  was  Dr.  Sung  H.  Ko 
( Code  2133 ) . 


Reviewed  and  Approved:  18  Nay  1990 


La^y  Freeman 
Associate/Technical  Director 
Surface  ASW  Directorate 


REPORT  DOCUMENTATION  PAGE 

form  Approved 

0MB  No.  0704-0188 

rteonif^  eura«n  lor  th'i  collection  of  mfornwiion  i>  ntimotM  10  eve'ige  <  nooc  Dcr  'twonM.  including  tn*  timo  for  roviowing  inttruciiom,  worcMna  oiittino  Oau  tourctl. 
ootherina  tnO  moinuining  the  date  needed,  and  completino  and  reviewing  the  collection  of  infer  motion  Send  comment!  rtaarding  ihn  burden  titifflatc  or  any  other  ateect  of  Ifm 
c^leaion  of  information,  including  tuggcttioni  for  reouemo  this  burden,  to  drauiinoton  Headduarters  ServKe!.  Directorate  for  information  Oberatiora  and  Kepom.  121$  jeffenon 
OayitHiphway.  Suite  1200.  Arlington,  va  222a2-o}02.andtotheOfficeof  Wanaoement  andeudoet.eaoerworkMeductionerotectfOTOO-OltW.Mrathington.  DC  20S0J 

'l.  ACENCV  USE  ONLY  a«dve  b/tni}  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

18  May  1990 

4.  TITLE  AND  SUBTITLE 

EVALUATION  OF  INTEGRALS  AND  SUMS  INVOLVING 
[si  n(Mx)/si  nlx)]*^ 

S.  FUNDING  NUMBERS 

PR:  A75215 

PE:  61152N 

C.  AUTNOR(S) 

Albert  H.  Nuttall 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOOR£SS(ES) 

Naval  Underwater  Systems  Center  ^ 

New  London  Laboratory 

New  London,  CT  06320 

B.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

TR  8689 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS{ES) 

Office  of  the  Chief  of  Naval  Research 

Arlington,  VA  22217-5000 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  OtSTRIBUTlON/AVAIUSILirY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTMCT  (Maximum  200  words) 

'^The  response  of  equispaced  arrays,  either  linear,  planar,  or  volumetric, 
to  distributed  spatial  fields,  typically  encounters  integrals  which  involve 
the  kernel  si n(Mx)/si n(x)  or  its  square.  Since  this  kernel  oscillates  rather 
fast  with  X  for  large  M  and  does  not  decay  with  x,  numerical  integration  of 
such  functions  can  be  very  time  consuming.  By  resorting  to  Parseval's 
theorem,  such  integrals  can  be  significantly  simplified,  requiring  only 
the  Fourier  transform  of  the  complementary  part  of  the  integrand.  This 
procedure  is  investigated  and  applied  to  several  typical  examples;  programs  .  - 

for  the  examples  are  also  included.  y/ /ay/'f 

1..  SUUICT  n»MS 

'  Summat 
Sine  f 
Arr;iy 

ation^  Fourie 

ion;  Parsev 

jnction,  Fast  F 

I’pfjnnnip.  Eouisr 

‘r  transform, 
al's  theorem, 

‘ourier  transform^ 
aced  arrav-  ( C.  r)\~— 

15.  NUMBER  OF  PAGES 

38 

16.  PRICE  CODE 

17.  SECURITY  classification 

OF  REPORT 

UNCLASSIFIED 

IS.  security  classification 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRAa 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

SAR 

NSN  7S40-01 -280-5500  St*n<J»fd  forfr  298  2-89) 

bt  ANS.  sio  I 


*:/ 


TR  8689 


TABLE  OF  CONTENTS 

Page 

LIST  OF  TABLES  ii 

LIST  OF  SYMBOLS  iii 

INTRODUCTION  1 

V. 

GENERAL  APPROACH  3 

Case  1  4 

Case  2  5 

Case  3  6 

Case  4  7 

EXAMPLES  9 

Example  A  9 

Example  B  11 

Example  C  14 

Special  Cases  16 

APPLICATION  TO  SUMS  19 

Case  n  =  0  20 

Case  n  «  1  21 

Case  n  -  2  22 

Example  22 

Some  Related  Sums  23 

.  SUMMARY  25 

APPENDIX  A  -  USE  OF  FAST  FOURIER  TRANSFORM  27 

APPENDIX  B  -  PROGRAMS  FOR  (6),  (10),  AND  (13B)  31 

REFERENCES  33 


i 


TR  8689 


LIST  OF  TABLES 

Table 

B-1.  Program  for  (6) 

B-2.  Program  for  (10) 

B-3.  Program  for  (13B) 


Accession  For 

NTIS  GRA&I 
OTIC  TAB 
Unannounced 
Justification. 


By. - 

Distribution/ 


1 

□ 


Availability  Codes 
[Avail  and/or 
Dlst  I  Special 


TR  8689 


LIST  OF  SYMBOLS 


N 

g(  t) 
G(  w) 
V 

H(W) 


Y 

h^(t) 

prime 

h2(t) 

H2(w) 

^2 
♦  (m) 

sub  r 


la 


2a 


V 


lb 


2b 


Ic 


2c 


8^(t) 

e 

Sj,(M,k) 


integer,  number  of  elements  in  line  array 
arbitrary  function  of  t,  (1) 

Fourier  transform  of  g(t),  (1) 

integral  of  product  of  g(t)  h*(t),  (2) 

Fourier  transform  of  h(t),  (2) 

argument  of  sine  functions,  (3) 

ratio  of  sine  functions,  (3) 

sum  with  prime  denotes  every  other  term,  (3) 

Fourier  transform  of  hj^(t),  (4) 

general  integral  (5) 

square  of  ratio  of  sine  functions,  (7) 

Fourier  transform  of  h2(t),  (8) 

general  integral  (9) 

autocorrelation  of  weights,  (12B) 

real  part,  (13B) 

integral  (18) 

integral  (22) 

integral  (27) 

integral  (28) 

integral  (33) 

integral  (34) 

increment  in  t,  (41),(A-1) 
impulse  train,  (41) 
convolution,  (42) 
discrete  sine  ratio,  (52) 


iii 


TR  8689 


INT 

V'{p) 
g(na) 
G(  w) 

N 

gc(nA) 

g(nA) 


integer  part,  (55) 
general  integral  (63) 
autocorrelation  of  ^(m),  (61) 
samples  of  g(t),  (A-1) 
approximation  to  G(a)),  (A-1) 
size  of  FFT,  (A-6) 
collapsed  version  of  g(nA),  (A-7) 
phase  modulated  g(nA),  (A-11) 
collapsed  version  of  §(nA),  (A-12) 


TR  8689 


EVALUATION  OF  INTEGRALS  AND  SUMS 
INVOLVING  (sin(Mx)/sin(x) ]” 

INTRODUCTION 

The  response  of  an  equiweighted  equispaced  line  array  to  a 
distributed  field  involves  the  kernel  sin(Mx)/sin(x)  or  its 
square,  depending  on  whether  the  voltage  or  power  response, 
respectively,  is  of  interest  (1;2].  Numerical  evaluation  of  such 
integrals  can  be  very  time  consuming  for  two  reasons;  this  kernel 
oscillates  quickly  with  x  for  large  M,  and  it  does  not  decay  with 
X.  This  necessitates  fine  sampling  in  x  and  large  integration 
regions,  both  of  which  can  lead  to  a  significant  computational 
burden,  especially  for  two-dimensional  or  three-dimensional 
arrays.  The  object  of  this  report  is  to  give  an  alternative 
numerical  procedure  that  can  be  very  advantageous  in  some  cases, 
and,  in  fact,  leads  to  closed  forms  for  some  examples. 

The  procedure  is  also  applied  to  summations  involving  the 
same  kernel.  Its  utility  depends  on  the  rate  of  decay  of  the 
complementary  part  of  the  original  integrand,  as  compared  with 
the  Fourier  transform  of  this  component.  In  any  event,  an 
alternative  is  presented  for  the  user  to  consider  in  any 
numerical  investigation. 


1/2 
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GENERAL  APPROACH 

For  arbitrary  function  g(t),  define  its  Fourier  transform  as 

G{(ti)  »  J  dt  exp(-i(»)t)  g(t)  .  (1) 

(Integrals  without  limits  are  over  the  range  of  nonzero 
integrand.)  Then  Parseval's  theorem  states  that  the  following 
two  alternative  integrals  are  equal: 

V  -  J  dt  g(t)  h*(t)  »  ^  J  dw  G(a>)  H*(a))  .  (2) 

Here,  HCw)  is  the  Fourier  transform  of  h(t).  Now,  if  H((»))  takes 
on  a  noticeably  simpler  form  than  h(t),  then  the  second  integral 
in  (2)  can  offer  an  attractive  alternative  to  the  first  integral 
in  (2).  That  will  indeed  be  the  case  here. 
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CASE  1 


For  integer  N  >  1  and  constant  r  >  0,  consider  the  special 
choice  of  h(t)  as 

M-1 

n«0 


M-1 

y~^  exp(iYtm)  , 


m=l-M 


(3) 


where  the  prime  on  the  latter  sum  denotes  skipping  every  other 
term.  Then  the  Fourier  transform,  according  to  (1),  is 


M-1 


m=l-M 


rm)  . 


(4) 


Substitution  of  (3)  and  (4)  in  (2)  yields 


V 


1 


dt  g(t) 


sin(MYt ) 
sin( Yt ) 


M-1 


E' 

m>=l-M 


G  (  Ym )  . 


(5) 


This  result  indicates  that  if  G(<i)),  the  Fourier  transform  of 
g(t),  can  be  evaluated,  then  the  t  integral  in  (5)  is  given  by  a 
finite  sum  of  equispaced  samples  of  G(w)  at  increment  2y.  The 
(complex)  function  g(t)  in  (5)  is  arbitrary,  except  that  the 
integral  must  converge.  When  GCoo)  cannot  be  analytically 
evaluated,  then  proper  application  of  a  fast  Fourier  transform 
procedure  to  g(t)  can  be  tailored  to  yield  precisely  the  equi¬ 
spaced  samples  required  for  the  right-hand  side  of  (5);  this 
technique  is  detailed  in  appendix  A. 
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If  function  Gita)  is  even  in  <0,  then  (5)  simplifies  to 


M-1 

2  ^  '  G(Yro)  for  M 
m-1 


2,4,6, 


M-1 

G(0)  +  2  ^ '  G(Ym)  for  M 
m-2 


1,3,5, 


(6) 


A  program  for  (6)  is  given  in  appendix  B, 


CASE  2 

For  integer  M  >  1  and  constant  y  >  0,  consider  the 
alternative  special  choice  of  h(t)  as 


M-1 


hjit) 


[flKT^]  -  E  -  1^'') 


n,  k-0 


M-1 


\  (M  -  |m|)  exp(i2Ytm)  , 
m-l-M 


(7B) 


where  we  used  (3).  There  is  no  prime  on  the  latter  summation 
because  all  terms  from  1-M  to  M-1  are  to  be  included.  The 
Fourier  transform  of  h2(t)  is 


H2(«) 


M-1 


2n 

m-l-M 


(M-  |m|)  6(w-2Ym)  . 


(8) 


The  use  of  (7A)  and  (8)  in  (2)  yields 
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^2  “ 


M-1 


9(t)  ^  (H  -  |m|)  G(2rm)  . 

m-l-M 


(9) 


Again,  the  integral  of  interest  is  given  by  a  finite  sum  of 
samples  of  the  Fourier  transform  of  g(t),  also  at  increment  2y 
in  «.  The  fast  Fourier  transform  technique  discussed  in 
appendix  A  is  relevant  here  also.  If  G(a>)  is  even  in  «,  then  we 
can  express  (9)  as 


M-1 


V2  “  M  G(0)  +  2 


y~'  (M  -  m)  G(2Ym)  for  all  M  >  1  .  (10) 

m=l 


A  program  for  (10)  is  given  in  appendix  B, 


CASE  3 


For  arbitrary  weights  and  frequencies  with 


exp(iYjj,t)  , 

m 


(llA) 


then  we  have  a  generalization  of  (3),  with 


83(00) 


w  &(oo  -  Y„)  . 
m  m 


(IIB) 


(Summations  without  limits  are  over  the  range  of  nonzero 
summand.)  Use  of  these  expressions  In  general  result  (2)  yields 


V 


3 


dt  g( t) 


★ 

w  exp(-ir„t) 
m  m 


G(y„)  . 
m 


(IIC) 


Again,  the  Fourier  transform  of  g(t)  is  required,  but  now  at 
general  arguments  {Yjjj}- 
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CASE  4 

Function  h2(t)  in  (7)  is  a  special  case  of  the  weighted  array 
power  response 

Wj^  exp(-i2Ytk)  «  ^  '  ♦(m)  exp(-i2Ytin)  ,  (12A) 

m 

where  ♦(m)  is  the  autocorrelation  of  the  weights; 


upon  use  of  (12A),  where  g(t)  can  be  complex  and  nonsymmetr ic . 
Thus,  integral  requires  the  autocorrelation  of  weights  {Wj^} 
and  the  Fourier  transform  of  g{t)  for  its  evaluation.  The 
earlier  result  in  (9)  corresponds  to  weights  “  1  for 
1  <  k  <  M. 

When  function  g(t)  is  real  (but  possibly  nonsymmetr i c )  and 
the  weights  are  real,  (13A)  can  be  simplified  to 
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V 


4 


♦  (0)  Gj.(0)  +  2 


E 


♦(m)  G^(2Ym)  , 


(13B) 


where  G^{oi)  is  the  real  part  of  Fourier  transform  G(w)  in  (1). 
A  program  for  (13B)  is  given  in  appendix  B. 
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EXAMPLES 


EXAMPLE  A 


The  first  example  of  interest  is 


g^Ct) 


(t-//)‘  +  8' 


,  P  >  0 


(14) 


Its  Fourier  transform  is 


Ga(w)  »  p  exp(-i/ya)-81  w|  )  , 


(15) 


for  which  the  real  part  is 


Gar(<»>)  “  p  cosifjui)  exp(-Plci)| 


(16) 


Since  integral  (5)  is  obviously  real  for  example  (14),  we 
obtain 


_ dt 

.  (t-|y)^  + 


_  sin(MYt) 
,2  sin(Yt) 


n— a 

E’ 

m«l-M 


Gar(Ym)  . 


(17) 


Substitution  of  (16)  in  (17)  yields  the  closed  form  result 


_ dt 

.  (t-^)^  + 


_  sin( MYt ) 
,2  sin(Yt) 


2n  r 

8D  M+3  ^M-1 


®M+1  ^M+1 


Ej^  (1  -  E2)  for  M  even 


2  (1  -  E^)  for  M  odd 


,  (18) 


where 


»  exp(-eYm)  ,  -  cos(^Ym)  ,  D  -  1  -  2  E2  C2  +  E^  .  (19) 
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A  program  for  (18)  and  (19)  follows;  it  is  written  in  BASIC  for 
the  Hewlett  Pacliard  9000  computer. 


10  INPUT  M, Bet a, Gamma, Mu  I  Beta  >  0,  Gamma  >  0 

20  B*Beta*Gamma 
30  C»Mu*Gamma 
40  E«EXP(-B*2) 

50  IP  (M  MODULO  2)-l  THEN  80 
60  F-COS(C)*SQR(E)*( 1-E) 

70  GOTO  90 
80  F-.5-.5*E*E 

90  A-E*COS(C*(M-l) )-COS(C*(M+l) ) 

100  A-A*EXP(-B*(M+1 ) )+F 

110  Vla«A*2*PI/(Beta*( l-2*E*COS(C*2 )+E*E) ) 

120  PRINT  Via 
130  END 


M-1 


E 

m«l-M 


(M 


|m()  cos(2/[/Ym)  exp( -2Py I m |  )  . 


(20) 


This  finite  sum  can  be  written  in  compact  form  by  use  of 
(3;  0.113],  Namely,  define  here 


E  -  exp(-2PY)  ,  C  =  cos(2//y)  ,  S  -  sin(2//Y)  , 

Epj  -  exp(-2BYM),  Cj^  =  cos(2/>yM),  =  sin(2/yYM), 

A-1-E^,  B=1+E^,  D*B-2EC.  (21) 

Then 
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2a 


dt 


(t-x/)‘  +  8' 


rsin(MYt)') 
Lsin(Yt)  J 


-  M[l -  K""  ®  ^  ^  =m]]  • 


A  program  for  (21)  and  (22)  follows. 


10  INPUT  N, Beta, Gamma /Mu  1  Beta  >  0,  Gamma  >  0 

20  Tb»2*Beta*Gamma 
30  Tm«2*Mu*Gamma 
40  E-EXP (-Tb) 

50  A»E*E 
60  B-l+A 
70  A-l-A 
80  C-COS(Tm) 

90  D-B-2*E*C 

100  Em-EXP(-Tb*M) 

110  T-(C*B-2*E)*(l-Em*COS(Tm*M) ) 

120  T-T+SIN(Tm)*A*Em*SIN(Tm*M) 

130  T-.5*M*A-T*E/D 

140  v2a-T*2*PI/(Beta*D) 

150  PRINT  V2a 
160  END 


EXAMPLE  B 

The  next  example  to  be  considered  is 

“  - 5 - 0  '  8  >  0  ,  a  >  0  .  (23) 

Since  9jj(t)  is  a  product  of  two  functions,  its  Fourier  transform 
Gj^(co)  is  given  by  a  convolution  of  the  individual  transforms. 

The  Fourier  transform  of  the  first  term  in  (23)  has  already  been 
encountered  in  (15),  and  the  Fourier  transform  of  the  second  term 
in  (23)  is  a  rectangle  located  on  interval  (-a,  a)  in  w. 
Therefore,  Gj^(w)  is  given  by  convolution 
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w+a 


Gj^(  w) 


2^  J  du  exp(-i/wu-e|ul  )  . 


(24) 


(O—a. 


Since  g^(t)  in  (23)  is  real,  we  need  only  evaluate  the  real  part, 

Sr 


(  w)  , 

of 

(w) 

.  With  the  aid  of 

auxiliary 

variables 

C  - 

fa) 

COS  (  //fa))  , 

s 

-  sinifjoi) 

'  s  - 

cosh(  Pfa)) , 

S  - 

—fa) 

sinh( 

Pfa)) , 

c  - 

a 

COS{  fJ 

a)  f 

s 

a 

•»  sin(//a) 

cosh( Pa) , 

s  * 

—a 

sinh( 

Pa)  f 

S  - 

c  c 

—fa)  fa) 

(P 

C 

a 

-  "  s.l  ♦ 

s  s 

—fa)  fa) 

(//  +  p 

S)' 

S  = 

s  c 

—a  a 

(P 

^fa) 

-  S„)  . 

c  s 

—a  a 

S)' 

find 

that  1 

Sr 

(  «) 

is  given 

by 

(25) 


Gbr(w)  = 


ap 


'P  -  exp(-Pa)  for  0  <  0)  <  a'l 

n 

(p  )  I  exp(-Pco)  B2  for  a  <  w 


(26) 


To  complete  the  description,  we  observe  that  Gbj.(«)  is  even  in  « 
because  gj2(t)  is  real.  A  program  for  Gbj.(w)  follows,  where  we 
have  made  the  following  identifications:  W  *  w,  A  ■  a,  B  ■  p, 

U  *  /y . 


10 

DEF  FNGbr(W,A,B,U) 

100 

20 

Wa-ABS(W) 

110 

30 

F-PI/(A*B*(B*B+U*U) ) 

120 

40 

Ea-EXP(-B*A) 

130 

50 

Ew-EXP(-B*Wa) 

140 

60 

Ca-COS(U*A) 

150 

70 

Cw-COS(U*Wa) 

160 

80 

Sa-SIN(U*A) 

170 

90 

Sw-SlN(U*Wa) 

180 

190 

IF  Wa<A  THEN  150 
Ra*l ./Ea 

T-(Ra-Ea)*Ca*(B*Cw-U*Sw) 

B2». 5*(T+(Ra+Ea)*Sa*(U*Cw+B*Sw) ) 
RETURN  F*Ew*B2 
Rw»l . /Ew 

T» ( Rw+Ew) *Cw* ( B*Ca-U*Sa ) 

Bl*. 5* (T+( Rw-Ew) *Sw* (U*Ca+B*Sa) ) 
RETURN  F*(B-Ea*Bl) 

FNEND 
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If  we  now  employ  (23)  in  (5),  we  obtain 


sin( at ) 
at 


sindlyt ) 
sin( Yt) 


' 


(27) 


m-l-M 


where  Gjj^((«))  is  given  by  (25),  (26),  and  its  even  property. 

Since  there  is  a  break  in  the  analytic  form  for  Gj^^(u))  at  w  «  ±a, 
it  is  not  reasonable  to  perform  the  summation  in  (27)  in  closed 
form;  those  terms  in  (27)  for  Y|m|  <  a  utilize  the  upper  line  of 
(26),  while  those  for  Y|n»|  >  a  utilize  the  lower  line  of  (26). 
However,  since  G|^^(w)  is  even  in  w,  the  simplification  in  (6)  is 
applicable . 

Instead,  when  (23)  is  substituted  in  (9),  there  follows 


(t-^) 


sin( at ) 
at 


rsin(MYt)] 
Lsin(Yt)  J 


(M  -  |m| )  G^^{2rm) , 


(28) 


where  Gjj^((«))  is  given  by  (25)  and  (26).  Again,  the  break  in  form 
of  Gj^^(w)  at  (0  »  +a  precludes  a  closed  form  result  for  the 
summation  in  (28);  also,  the  simplification  in  (10)  is 
immediately  applicable  to  (28). 
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EXAMPLE  C 


The  final  example  is 

”  - 5 - 2  '  P>0'  o>0.  (29) 

The  Fourier  transform  of  the  second  term  in  (29)  is  a  triangle 
located  on  interval  (-2a,  2a)  in  w.  Therefore,  0^(0))  is  given  by 
convolution 

(«)'•- 2  a 

Gc(w)  ”  2^  J  exp(-iA/u-Plu|  )  ^1  -  .  (30) 

w-2a 


Because  g^(t)  is  real,  only  the  real  part  of  (30)  is  needed. 
This  tedious  calculation  has  been  carried  through,  with  the 
following  result;  define  auxiliary  variables 


R  » 

fl2  2  ^ 

8  -//  ,  I  = 

28/^,  D 

2 

^  fi  +fj  , 

exp(-8(*))  f 

C 

—a 

-  cosh( 28a) , 

S  «  sinh(28a), 

•“Ot 

C  - 
a 

cos  (  2/ja)  , 

3 

Ul 

-  cosh(8<i»)  , 

«  sinh( 8«)  / 

cos(^co)  , 

Cl 

-  C  C  (R  C 

-IS) 

+  s,  S  ^ 

(R  S 

+  I  C  )  , 

1 

a 

3 

3 

1 

a 

—CO  CO 

a 

a 

Co 

-  C  C  (R  C 

-IS) 

+  S  S 

(R  S„ 

+  I  c,  )  , 

2 

—a  a  w 

CO 

—a  a 

CO 

CO 

C  - 

R  C  -IS 

CO  CO 

• 

E^  «  exp( -28a) , 
S  >  sin(2;/a)  , 

(X 

-  sinifjoi)  , 


(31) 


Then  we  find  that 


I'D  8  (2a  -  w)  -  C  +  for  0  <  w  <  2a'j 


G.  .  (w) 


2a^8D' 


-  E  C  +  E  C-  for  2a  <  w 
0)  w  2 


(32) 
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Also,  G  (w)  is  even  in  w.  A  program  for  G  (w)  is  listed  below, 

Cl  Cl 

where  W  ■  w,  A  ■  a,  B  ■  |3,  U  ■  //. 


10  DBF  FNGcr(W,A,B,U)  1  A  >  0  ,  B  >  0 

20  Wa-ABS(W) 

30  Tb-2.*A*B 

40  Tu-2.*A*U 

50  Bw-B*Wa 
60  Uw-U*Wa 
70  B2-B*B 
80  U2-U*U 
90  R-B2-U2 

100  I-2.*B*U 

110  D-B2-I-U2 

120  Ew-EXP(-Bw) 

130  Ea-EXP(-Tb) 

140  Ca-COS(Tu) 

150  Sa-SIN(Tu) 

160  Cw-COS(Uw) 

170  Sw-SIN(Uw) 

180  C-R*Cw-I*Sw 

190  IF  Wa<2.*A  THEN  250 

200  Ra-l./Ea 

210  C2-.5*(Ra+Ea)*Ca*C 

220  C2-.5*(Ra-Ea)*Sa*(R*Sw+I*Cw)+C2 

230  T-Ew*(C2-C) 

240  GOTO  290 
250  Rw-I./Ew 

260  Cl-.5*(Rw+Ew)*Cw*(R*Ca-I*Sa) 

270  Cl-.5*(Rw-Ew)*Sw*(R*Sa+I*Ca)+Cl 
280  T-D*(Tb-Bw)-Ew*C+Ea*Cl 
290  RETURN  PI *T/( Tb*A*D*D ) 

300  FNEND 


We  now  substitute  (29)  into  (5)  and  get 


Ic 


dt 


rsin( at)]  sin(MYt ) 
L  at  J  sin(Yt) 


M-1 

E 

m-l-M 


n— X 


(Ym)  ,  (33) 


where  G^^(w)  is  given  by  (31),  (32),  and  its  even  character.  The 
brea)(  in  form  in  (32)  at  co  «  ±2a  precludes  a  closed  form  for  the 
sum  in  (33).  However,  (6)  is  still  applicable. 
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When  (29)  is  utilized  in  (9),  there  follows 


dt 

■sin(  at )  ■ 

2 

rsin(MYt )] 

(t-//)^  +  8^ 

L  at  J 

Lsin(Yt)  J 

M-1 

-  ^  (M  -  |m|  )  G^j.(2Ym)  . 

m-l-M 

Equation  (10)  may  also  be  employed  here. 


SPECIAL  CASES 


If  we  set  M  =  1  in  (17),  there  follows 


2^^.  «2  '  *^ar^°^  “  8  ' 


+  8 


where  we  used  (16).  The  same  case  in  (27)  yields 


dt 


+  8 


2  . 


{8  -  exp(-8a)  18  cos((/a)  -  /j  sin(A/a)]}  , 


upon  use  of  (26)  and  (25).  Finally,  from  (33), 


o8(8^+/t/^) 


(34) 


(35) 


(36) 
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_ dt _ 


^sin(  at) 


cr 


(0)  - 


- 5 -  f2ae(e2  +  //2)  -  R  +  (R  -  I  S^))  , 

2a^e(a^+//^) 

using  (32)  and  (31). 


(37) 
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APPLICATION  TO  SUMS 


In  this  section,  it  is  more  convenient  to  use  Parseval's 
theorem  (2)  in  the  form 

V  -  J  dt  g(t)  h*(t)  -  J  df  G(f)  H*(f)  , 

where  Fourier  transform 


(38) 


G(f)  -  J  dt  exp(-i2iift)  g(t)  . 


(39) 


Now,  we  take  as  our  candidate  h(t)  function, 


h(t)  =  p(t)  AS*(t)  , 


(40) 


where  S*(t)  is  the  infinite  impulse  train 


S^(t) 


E 


8(t  -  kA)  . 


(41) 


The  Fourier  transform  of  h(t)  is  then 


H(£)  .  P(£)  e  8j/^(£)  .  ^  p(£  -  I)  , 


(42) 


where  P(f)  is  the  Fourier  transform  of  p(t),  9  denotes 
convolution,  and  we  have  utilized  the  fact  that  the  Fourier 
transform  of  impulse  train  A5^(t)  is  another  impulse  train, 

Substitution  of  (40)  and  (42)  in  (38)  yields 
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E 


g(kA)  p  (kA) 


EJ 


df  G(f)  P 


(43) 


For  general  p(t)  and  P(f),  this  will  not  be  a  useful  relation, 
since  the  right-hand  side  of  (43)  is  an  infinite  sum  of 
integrals.  However,  we  will  be  interested  here  only  in  the 
special  cases  of 


P**"’  ■  [sln(rt)’]  '  "  • 


(44) 


CASE  n  =  0 

For  n  =  0,  the  above  relations  specialize  to 

p(t)  =  1  ,  P(f)  =  S(f)  , 


H(f)  = 

vo-E 

k 


(45) 


This  is  a  discrete  version  of  Parseval's  theorem.  Although  one 
infinite  sum  has  been  traded  for  another,  we  can  now  choose  that 
alternative  that  has  the  most  rapidly  decaying  (and/or  easily 
computed)  summand  for  numerical  evaluation. 
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(46) 


(47) 


(48) 


Again,  we  have  an  alternative  infinite  sum  (48)  that 
hopefully  decays  faster  than  the  original  sum  (47).  The 
sin( Mx )/sin( X )  term  does  not  help  convergence  in  (47)  because 
this  term  never  decays  for  large  x.  Although  (48)  is  a  double 
sum,  the  summation  on  m  only  contains  M  terms;  the  utility  of 
(48)  depends  heavily  on  the  asymptotic  decay  of  G(f)  for  large  f. 
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CASE  n  »  2 

With  the  aid  of  (7),  we  now  find 


M-1 


“  [fl^Ttr]  “  ^  (M  -  |m|)  exp(i2Ytm)  , 

m-l-M 


M-1 


P{f)  =  >  (M 

M-1 


(M  -  |m|  )  &  ^f  -  , 


H(f)  =  ^  X]  i  “  ^]  ' 


k  m=l-M 


k 

M-1 


E  E  • 


k  in=l-M 


(49) 


(50) 


EXAMPLE 


Consider,  as  in  (14)  and  (16), 


g^(t) 


2  2  ' 


“  fl  cos(2ii^f)  exp( -2ii6  I  f  I  ) 

or  p 


(51) 


The  summations  in  (47)  and  (49)  are  very  slowly  decaying,  leading 
to  difficulty  in  attaining  accurate  results.  The  alternatives  in 
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(48)  and  (50),  on  the  other  hand,  have  exponential  decay  and  can 
be  evaluated  quite  accurately.  The  additional  examples  given 
earlier  in  (23)-(26)  and  in  (29)-(32),  along  with  the 
corresponding  programs,  lend  reasonable  alternatives  to  some 
otherwise  lengthy  numerical  calculations. 
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The  sum  over  a  double  interval  is 


2N-1 

rO  for  M  -  0,2,4, . 

y  s„(«,k)  .  . 

o 

1 

UN  for  M  =  1,3,  ..  .  ,2N-lJ 

The  correlation  on  the  second  variable  of  S.,  is 

N 

N-l 

y  Sj^(M,k)  Sj^(M,k+j)  =  N  Sjj(M,j)  for  0  <  M  <  N  and 

k=0 

Finally,  the  correlation  on  the  first  variable  is 


N-1 

Sj^(M,k)  Sj^(M+2L,k) 
'M(N  -  M  -  2L) 

+  . 

N(3M  +  2L  -2N)  -  M(M  +  2L) 

for  all  M,  L,  N,  where 


M(M  +  2L)  + 


for  0  <  M  +  L  <  N' 
for  N  <  M  +  L 


M  =  M  MOD  N  ,  L  =  L  MOD  N  . 


(56) 


.  (57) 


(58) 


(59) 
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SUMMARY 


Extensions  to  integrals  involving  [ sin(Mx)/sin(x) ]”  for  n  >  2 
are  possible,  based  upon  the  results  presented  here.  For 
example,  starting  from  (12A)  for  arbitrary  weights,  we  could 
consider 


where 


hjCt) 


h^t) 


ij;(p)  exp(-i2Ytp)  , 


(60) 


V'(p) 


♦  (m) 


♦  (m-p) 


(61) 


is  the  autocorrelation  of  sequence  {♦(m)}  defined  in  (12B). 
Therefore,  Fourier  transform 


H5(w) 


8(co  +  2yP)  t 


giving  rise  to 


(62) 


Vg  =  J  dt  g(t)  hgCt)  =  ^  V/(p)  G(2yP)  •  (63) 

P 

The  case  of  equal  weights  in  (12A)  now  corresponds  to  n  »  4 

in  the  sine  function  ratio  above,  and  ^Kp)  is  the  autocorrelation 
of  a  triangular  sequence. 
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The  evaluation  of  integrals  and  sums  involving  the  term 
I sin(Mx)/sin(x) can  often  be  simplified  by  the  use  of 
Parseval's  theorem  because  this  term  has  a  Fourier  transform 
which  is  a  finite  sum  of  delta  functions.  Major  effort  can  then 
be  concentrated  on  getting  the  Fourier  transform  of  the 
complementary  part  of  the  integrand.  This  procedure  has  been 
applied  here  to  several  examples  which  arise  in  evaluation  of  the 
response  of  equispaced  arrays  to  distributed  spatial  fields.  For 
more  complicated  fields,  a  fast  Fourier  transform  procedure 
combined  with  the  above  result  leads  to  a  very  efficient  method 
of  integral  evaluations.  Applications  of  this  procedure  have 
been  made  in  [ 5 ] . 
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APPENDIX  A  -  USE  OF  FAST  FOURIER  TRANSFORM 

The  summations  for  and  V2  in  (5)  and  (9),  respectively, 
require  the  evaluation  of  the  Fourier  transform  of  g(t),  namely 
G(oi>),  at  equispaced  increment  2y.  But  this  latter  function  can 
be  approximated  by  means  of  the  trapezoidal  rule  according  to 

G(w)  «  J  dt  exp(-iMt)  g(t)  * 

=  L  ^  ‘  exp(-i(«)An)  g(nA)  ■  G(w)  «  ^  G^w  -  ,  (A-1) 

n  n 

where  A  is  the  sampling  increment  in  t.  The  latter  summation  in 
(A-1)  indicates  aliasing  lobes  separated  by  2n/t  on  the  co  axis. 

In  order  to  control  aliasing,  we  must  choose  A  small  enough,  say 
A  <  Aq.  Then  samples  of  approximation  G(«)  in  (A-1)  at  multiples 
of  2y  are  given  by 


G(  2 Yin)  »=  A 


exp(-i2YAmn)  g(nA)  . 


{Pi-2) 


Now  since  A  is  arbitrary,  except  for  upper  limit  A^,  choose 


Ny  ' 


(A-3) 


where  N  is  an  integer  and  y  is  the  prescribed  increment  in  a>.  In 
order  that  A  be  less  than  A^,  we  must  ta)te  integer 


^"o  ‘ 


(A-4) 


Use  of  (A-3)  in  (A-2)  gives 
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G{2rm)  -  A 


exp(-i2iimn/N)  g(nA) 


(A-5) 


X 

A  ^  '  exp(-i2nmn/N)  g^(nA)  , 


(A-6) 


where  "collapsed"  sequence  [4;  pages  4-5] 


g^(nA)  ■  y  '  g(nA  +  kNA)  for  0  <  n  <  N  -  1  . 


(A-7) 


The  manipulation  from  (A-5)  to  (A-6)  is  exact;  it  avoids 
truncation  error  normally  associated  with  functions  g(t)  which 
decay  slowly  with  t.  The  sum  on  k  in  (A-7)  must  be  carried  out 
(for  each  n)  until  negligible  values  for  g  are  encountered  for 
both  positive  as  well  as  negative  values  of  k. 

Equation  (A-6)  indicates  that  values  of  G(2Yni)  for  m  »  0  to 
N  -  1  are  available  by  an  N-point  fast  Fourier  transform  when  N 
is  a  power  of  2.  Values  for  negative  m  are  available  in  location 
m  mod  N.  In  order  to  get  all  the  desired  values  of  G(2Ym) 
required  for  (9),  without  aliasing,  we  also  require  that  N/2  >  M. 
Thus,  the  final  condition  on  integer  N  is 


N  >  max^^^,  2Mj  . 


(A-8) 


For  the  case  of  (5),  where  the  increment  on  w  is  2y,  but 
starting  at  to  *  Yf  we  return  to  (A-1)  to  find  that 


G  (  Y  +  2  Yin ) 


exp(-iY6n  -  i2YAmn)  g(nA) 


(A-9) 
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The  same  choice  of  A  in  (A-3)  now  yields 


G(  Y  +  Zym)  “  A  )  exp(-i2iimn/N)  exp(-iiin/N)  g(nA)  .  (A-10) 


This  result  is  identical  to  (A-5)  except  that  g(nA)  must  be 
replaced  by 


exp(-iiin/N)  g(nA)  ■  g(nA)  . 


(A-11) 


Calculation  of  the  collapsed  version  of  g  is  eased  by  the 
observation  that 


g^(nA)  s  9(nA  +  kNA) 


exp(-in(n  +  kN)/N)  g(nA  +  kNA)  >= 


exp(-inn/N) 


1)  g(nA  +  kNA)  for  0  <  n  <  N  -  1  ,  (A-12) 


thereby  leading  exactly  to 


G  (  Y  +  2  yhi  )  *  A 


y  '  exp( -i2!imn/N)  g^(nA) 


(A-13) 


The  leading  phase  factor  in  (A-12)  only  needs  to  be  evaluated  at 
N  different  values,  and  the  sum  in  (A-12)  requires  differencing 
of  "adjacent"  samples  of  g  spaced  by  NA,  rather  than  the  straight 
summation  previously  adequate  for  (A-6)  and  (A-7).  Condition 
(A-8)  applies  here  as  well. 
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APPENDIX  B  -  PROGRAMS  FOR  (6),  (10),  AND  (13B) 


Table  B-1.  Program  for  (6) 

10  M-7  1  >  0 

20  Gamma*1.31  !  >  0 

30  DOUBLE  M,MS  !  INTEGERS 

40  S-0. 

50  IF  (M  MODULO  2)-l  THEN  110 
60  FOR  Ms-1  TO  M-1  STEP  2 
70  S-S+FNG(Gamma*Ms) 

80  NEXT  Ms 

90  Vl-2.*S 

100  GOTO  150 

110  FOR  Ms-2  TO  M-1  STEP  2 
120  S-S+FNG(Gamma*Ms ) 

130  NEXT  Ms 

140  Vl-FNG( 0. )+2. *S 

150  PRINT  M, Gamma, VI 

160  END 

170  I 

180  DEF  FNG(W) 


Table  B-2.  Program  for 

(10) 

10 

M-6  ! 

>  0 

20 

Gamma-. 71  1 

>  0 

30 

DOUBLE  M,NS  ! 

INTEGERS 

40 

G2-2 . * Gamma 

50 

S-0. 

60 

FOR  Ms-1  TO  M-1 

70 

S-S+(M-Ms ) *FNG(G2*MS) 

80 

NEXT  Ms 

90 

V2«M*FNG(0. )+2.*S 

100 

PRINT  M, Gamma, V2 

110 

END 

120 

1 

130 

DEF  FNG(W) 
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Table  B-3.  Program  for  (13B) 


10 

M-9 

! 

>  0 

20 

Gamma- .79 

t 

• 

>  0 

30 

DOUBLE  N,Ms,Ks 

f 

• 

INTEGERS 

40 

DIM  W(IOO) 

50 

REDIM  W(1:M) 

60 

CALL  Weights(M,W( * ) ) 

! 

REAL  WEIGHTS 

70 

G2-2 . *Gamma 

80 

S-0. 

90 

FOR  Ms-1  TO  M-1 

100 

Phi-0. 

110 

FOR  Ks-Ms+1  TO  M 

120 

Phi-Phi+W( Ks ) *W( Ks-Ms ) 

t 

CORRELATION  OF  WEIGHTS 

130 

NEXT  Ks 

140 

S-S+Phi*FNGr(G2*Ms) 

150 

NEXT  Ms 

160 

Phi-0 . 

170 

FOR  Ks-1  TO  M 

180 

Phi-Phi+W( Ks ) *W( Ks ) 

190 

NEXT  Ks 

200 

V4-Phi*FNGr(0. )+2.*S 

210 

PRINT  M, Gamma, V4 

220 

END 

230 

! 

240 

SUB  Weights (DOUBLE  M,REAL 

W(*)  ) 

250 

DOUBLE  Ks 

1 

INTEGER 

260 

T-2. *PI/M 

270 

FOR  Ks-1  TO  M 

280 

D-KS-.5 

290 

W(Ks)-l. 

{ 

FLAT  WEIGHTS 

300 

W(Ks)-.5-.5*COS(T*D) 

1 

HANN  WEIGHTS 

310 

W(Ks)-.54-.46*COS(T*D) 

1 

• 

HAMMING  WEIGHTS 

320 

NEXT  Ks 

330 

MAT  W-W/SUM(W) 

1 

NORMALIZATION 

340 

SUBEND 

350 

1 

360 

DEF  FNGr(W) 
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